#pragma once

void print_matrix(double** phi, int n, char* buf);

void print_matrix2(double** phi, double** eta, int n, char* buf);

void print_matrix3(double** phi, double** eta, double **xx, int n, char* buf);

void print_vector(double *x, int n, char* buf);

void Householder(double* x, double *v, double *beta, double n);

void HouseholderQR(double **A, double Q, int m, int n);

void matrix_product(double **a, int row_la, int row_ra, int col_la, int col_ra,
	                double **b, int row_lb, int row_rb, int col_lb, int col_rb,
	                double **c, int row_lc, int row_rc, int col_lc, int col_rc)

